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Abstract 

This note constitutes a corrie endum to the article of lAzomahoul (|2009l '). JSPI, 139:2581- 
2597. The aggregation of isotropic four nearest neighbors autoregressive models on the lattice 
1? , with random coefficient, is investigated. The spectral density of the resulting random 
field is studied in details for a large class of law of the AR coefficient. Depending on this law, 
the aggregated field may exhibit short memory or isotropic long memory. 
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1 Introduction 



Contemporaneous aggregation of random parameter autoregressive processes is a well - known 
scheme to produce long memory in time seri e s. Th is idea first appeared in Robinsoiil (1978) 
and was more deeply investig ated in iGrangeii (If 9801) . A lot of papers have then been d evote d 
to this subject, among them Gonpalvez and Gourieroux ( f988[ ). Qppenheim and Viano ( 2004 ). 
Zaffaronil ( 20041 ). This construction of long memory is famous in econometrics, since it illustrated 



the fact that long memory may arise in macroeconomic time series, although micro-dynamics are 
short range dependent. 

More recently, the concept of long memory has been extended to random fields on a lattice. 
In this setti ng, the intensity of long memory may depend on the direction or may be isotropic. 
We refer to iLavancie 3 (l2005l) for a revi e w. From a statisti cal p oint of view, the isotropic case is 
investigated in Anh and Lunnev| (1995 ). Frias et al. (2008) and Wane (l2009l). Som e anisotropic 



models are studied in lBoissv et al.l (|2005l) . Guo et al. (200£ ) and Beran et alT ( 2009[ ). A statistical 
test to detect long memory in random fields is provided in Lavancieil (j2008h ' 

As for time series, it seems natural to study the effect of aggregation in random fields, in 
particular the emergence of long memory. The aggregation in this case in not contem 
but is implemented at each point of the lattice. This scheme has been described in 
2005 ). where it is applied to produce some anisotropic models of long memory. In lAzomahou 



poraneous 



Lavanciei 



2OO90, aggregation of isotropic four nearest neighbors autoregressive models on the lattice Z is 
studied. These autoreg ressive models are very natural models on 1? and their introduction dates 
back to Whittle! ( 1954). Their aggregati on is expected to produce isotropic long memory in some 
cases, as announced in Azomahou ( 2009f ). Unfortunately, this article contains false results and the 
emergence of isotropic long memory by aggregation is not clearly descr ibed. 

The present note is a corrigendum to the article of lAzomahorJ ( 20091) . In section [51 the general 
scheme of aggregation for random fields on 1? is recalled. In section[3J the aggregation of isotropic 
four nearest neighbors autoregressive models is considered. First, we recall the necessary and 
sufficient conditi ons for existen ce of a sta tionary s olutio n to such autoregressive equations, which 
are known since IWhittld (|l954h (see e.g. iGuvonI (Il995h l. but are improperly stated in Lemma 1 
in Azomahoul (j2009l) . Then, the existence of a stationary aggregated field in is studied. Under 



some general conditions on the law of the autoregressive parameter, this existence is ensured in 
Theorem [T] The behavior of the spectral density of the induced field is also presented. The 
memory properties of the aggregated f i eld ar e then deduced. These results are a corrected version 
of Propositions 2 and 3 in Azomahoul ( 2009f ). 



2 Aggregating autoregressive random fields with random co- 
efficients 

The definition o f an autoregressive fie l d in 1? is very similar t o the definition of an autoregressive 
time series (see Brockwell and David (Il991)) and appears in Whittld (1954). Basic results are 
presented in iRosenblattI (|l985[ ) and lGuvonI (|l995f ). 

Let L\ (resp. Lq) be the lag operator acting on the first (resp. second) index of the field 



i.e. LiXij — Xi. 



and L2Xi 



X, 



The field {Xi_j)(^i j-)^j2 is an autoregressive field if there exists a white noise {£i,j)(ij)^z^ in 
and a complex function P{zi,Z2) = '^(^k i)i£S '^'^''■^1^2^ where S* is a finite subset of 1? and 
(flfc,i)(fc,i)GS ^6^1 coefficients, such that for all G 



P{Li, L2)Xij — £, 



(2.1) 



It is well known th at such a n equa ti on admits a statio nar y solution if an d only if P(e*^% e*^^) 7^ 



for all (6*1, 6*2) (see IWhittl^ (|l954f) . IRosenblattI (|l985l ) or iGuvonI (|l995l )). In this case, the field 
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{Xij)(ij-jf=z^ may be written 

Xij = ^ Ck,iSi-k,j-i (2.2) 

where (cfe,;) are the coefficients involved in the Laurent expansion of P~^{zi, Z2) in the neighbor- 
hood of {|zi| — 1, 1^2! = 1| and satisfy \ck i\ < 00. Note that the series in (|2.2p converges almost 



surely and in (see Brockwell and David ( 199lh . Proposition 3.1.1, for a proof in the time series 



setting which can be easily adapted to random fields). 

Let us now turn to autoregressive fields with random coefhcients. They are defined as in 
(|2.ip where the coefficients {o,k,i){k,i)£S defining P are random. Assuming that for almost every 
(afcj)j -P(e'^Se*^^) ^ for all (01,62), we deduce easily that the representation (12.21) (whe re the 



coefficients {ck,i) are now random) holds almost surely. See lOppenheim and Viand (|2004l ) for a 
justification in the time series setting which remains valid for random fields. However, for (|2.2p to 
be true in L^, it is necessary and sufficient that 



E [ \P{e'^\e'^^)\-^dXidX2 < 00, 

J[-7r,7rl2 



(2.3) 



where the expectation is with respect to the law of the coefficients (afc.i)(fc,i)GS (^^^ the same 
references). If (|2.3p holds, the spectral density of {Xij)^i jj^i2 is 

fx{XuX2)^^E\P{e^^\e^^-^)\-\ (2.4) 

where ct^ denotes the variance of the white noise (ci j")(i j)gz2 . 

Finally, aggregation of random fields is implemented according to the classical scheme for time 
series. Consider N independent copies of autoregressive fields with random coefficients as above, 
i.e. for 71 = 1, . . . , A^, for all (i, j) e Z^, 

p(")(Li,L2)4")=eg, (2.5) 

where (e^")) are i.i.d. copies of white noises (called idiosyncratic components) and the p(")'s come 
from i.i.d. copies of the random coefficients iak,i){k,i)es- Aggregation of these A^ fields defines the 
field Z(^) as, for aU £ 1? , 

N 



^E^ff- (2-6) 



^ 1 

n— 1 



Under condition (12.31) . since the (Xj'-"'')„'s are i.i.d. the field Z^^^ shares the same 1? properties 
as the (x|"'')„'s. In particular, its spectral density is (|2.4p . When A^ goes to infinity, from the 
central limit theorem, Z'--^^ tends to a Gaussian field Z. The aggregated field Z is thus a zero 
mean Gaussian field with spectral density (|2.4p . Moreover, contrary to Z'^^\ Z is ergodic (see 
R obinsonI (119781) 1. 



3 Isotropic long memory random fields as aggregated fields 

In this secti on we study the p articular aggregation of isotropic four nearest neighbors autoregressive 
fields as in lAzomahoul (|2009i ). Unfortunately that article contains false statements. Our results 
correct them and provide further properties. 

For each n G N, AT^"' is the autoregressive field following the representation, for all (i, j) G 1? , 



where e*^"^ is a white noise in L^. This representation corresponds to (|2.ip where PjL i, L2) = 
1 — 9{Li+L^^ +L2 + L2^)- When 9 is not random, this model has been first considered in I Whittle 
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(Il954l) and is stu died for instance in iGuvon (|l995[ l . Contrary to the statement of Lemma 1 in 



Azomahou ( 2009h . we deduce from these references (or by easy computation) that the condition 



\9\ < J is necessary and sufficient so that p.ip admits a stationary solution in 

Let us now assume that 6 is random and belongs almost surely to [0, j[ (see Remark [5] for the 
case 9 < 0). More precisely, let us suppose that the law of 9 admits a density of probability of the 
form (f>(a;)(i; ~ x)" , for < a: < ;| , where a > — 1 and $ is some function whose properties will 
be specified later. According to (|2.4p . if (|2.3p holds, the spectral density of X*^"^ defined by p.ip 
is in this case 

Jx(M,^2) - _ 22.(eos(Ai) + cos(A2)))2 

Assumption (|2.3p is equivalent to the integrability of fx on [— tt, tt]^. Conditions on a and $ 
which imply (j2.3p are deduced from the following lemma. 

Lemma 1. Let a > —1 and ^ be a bounded function on [0, j], continuous at j with $(^) ^ 0. 
Consider, for any cr^ > 0, the function f defined on [— 7r,7r]^ by 



Hx)a-xr 



/ : (Ai,A2) ^ / , , ;:\ '———-rdx. (3.2) 

V i, 47r2 7o (1 - 2a;(cos(Ai) +cos(A2)))2 ^ ^ 

('i^ w/ien a> 1, f is continuous at any point (Ai, A2) G [— tt, Tr]^, 

(ii) when —1 < ct <1, f is continuous at any (Ai,A2) 7^ (OjO) in [— 7r,7r]2 and diverges to +00 
when (Ai,A2) (0,0). Moreover, as (Ai,A2) -J' (0,0), 

- < a < 1, 

/(Ai,A2)^c„(A? + A^)"-i (3.3) 
c„ = ^16-"$(i) /+- ^^du. 

— if a = I and assuming ^ is a P- Holder function (j3 > 0), 

/(Ai,A2)^ci|ln(A? + A2)| (3.4) 

where ci = ^16-i$(i). 



The proof of this lemma is postponed to the appendix. As a consequence of Lemma [1] (|2.3|) 
holds if a > and $ satisfies the assumptions in Lemma [TJ In this case, the aggregated field 
Z, constructed as the limit of (|2.6p when N — > 00, belongs to and its spectral density is 
p.2p . The memory properties of Z are deduced from the behavior of its spectral density. When 
it is continuous, then Z exhibits short memory. When the spectral density is unbounded, then 
the autocovariance function of Z is not summable and Z is said to exhibit long memory. These 
properties are summarized in the foll owing theor e m. It provides the spectral properties improperly 
presented in Propositions 2 and 3 of Azomahou ( 2009f ). 



Theorem 1. Let X^'^\ n £ N, be a sequence of i.i.d isotropic autoregressive fields defined by 
p.ip . where 9 is random with probability density ^{x){j — x)" , < a; < -j. Assume that a > 
and that ^ is a bounded function on [0,j], continuous at j with ^{j) 7^ 0. If a = 1, assume 
moreover that ^ is a (3 -Holder function with /3 > 0. Let Z be the aggregated field obtained as the 
limit of (|2.6p when N ^ 00. Then Z is a zero mean Gaussian random field in L^ with spectral 
density (13. 2p . where denotes the variance of the white noises in p.ip . In particular, 

• if a > 1, Z is short-range dependent, for its spectral density is continuous everywhere, 

• ifO<a<l,Z exhibits long memory and its spectral density behaves as p.3p at 0, 

• if a — 1, Z exhibits long memory and its spectral density behaves as p.4p at 0. 
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According to the resuhs recalled in section [51 the proof of this theorem is a direct consequence 
of Lemma [TJ 

Remark 1. Following Definition 1 in Lavancie"^ ( 2005 ) . a random field over exhibits isotropic long 



memory if its spectral density behaves at as \X\'^L b ^]^^ j ~2 < 7 < 0, where A — (Ai, A2), 

|.| denotes the euclidean norm, L is slowly varying at infinity and where b is continuous on the 
unit sphere in M^. According to (|3.3p . this is the case for Z when < a < 1 with 7 = 2a — 2. 
When a = 1, the behavior of the spectral density at is also a diverging function of |A|. For this 
reason, one may also say that Z exhibits isotropic long memory, even if it does not suit the latter 
definition. 

Remark 2. From a general point of view, the law of 9 may be supported on the whole interval 
] — i , i [. From Theorem [I] the behavior of the spectral density at is driven by the behavior of 
the probability density of 9 near j. It is easy to see that under similar conditions, the behavior 
of the probability density of at —j rules the behavior of the spectral density at points (Ai, A2) 
such that |Ai| = IA2I = tt. In particular, some similar spectral sing ularities may arise at th ese 
points, leading to the so-called seasonal long memory phenomena (see Qppenheim et all (j2000f) in 
the time-series setting). 

A Proof of Lemma [T] 

For any fixed x < i, the integrand in (I3.2p is continous with respect to (Ai,A2) € [— tt, tt]^. 

Moreover, since cos(Ai) -I- cos(A2) < 2 , it is bounded by ^^{j — which is integrable when 

a > 1. So (i) is proved thanks to the dominated convergence theorem. 

Similarly, for any (Ai, A2) in a sufficiently small neighborhood of some fixed (Ai, A2) 7^ 0, there 
exists e > such that cos(Ai) + cos(A2) < 2 — 2e. Thus for any fixed x < ^, for any such (Ai, A2), 
the integrand in (j3.2l) is bounded by e~'^^{x){j — x)". This gives the continuity of / at any point 
(Ai, A2) 7^ when a > —1. 

On the other hand, changing j — x into z in p.2p leads to 



/(Ai,A2) 



47r2 (2 - cos(Ai) - cos(A2))2 Jo cos(Ai)+cos(A2) ^ ^ 

2— cos(Ai ) — cos(A2) 



Letting Ax — 2 - cos'^( a 1")"- cos'^( ) ^^^^ applying the change of variables u — 4:zAx, we obtain 

^^ ^ £! |cOs(Ai)+COs(A2)|-° f^^ Wmlil-f^)) 

^' 47r2 (2-cos(Ai)-cos(A2))i-" cos(Ai) + cos(A2) io (l + u)^ 

(A.l) 

Note that A\ — )■ -t-oo as (Ai, A2) — (0, 0). Since $ is continous at j and bounded on [0, |], the 
dominated convergence theorem applies when — 1 < a < 1 and easy computations give 
When a = 1, we get from (lA.ll) . 



/(Ai,A2)^-^ ——ih+h), (A.2) 

47r^ (cos(Ai) + cos(A2))^ 



where 



/r = r (-) /2 ^ r 

Assuming $ is a /3-Holder function on [0, |] when a = 1, we have 



4 I Ax 



\h\ < 



\4A^ 1/3 



AT Jo 



{1 + uy 



du, 



du. 



(A.3) 
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where c > 0, which shows that /2 — > as (Ai, A2) — > (0, 0). For Ii, note that A\ > if (Ai, A2) is 
small enough, so direct computation leads to 



The equivalence ((5^ then follows easily from (|A.2p . (|A.3p and (|A.4p . 
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